Random walk of a swimmer in low- Reynolds- number conditions 
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Swimming at a micrometer scale demands particular strategies. When inertia is negligible com- 
pared to viscous forces, hydrodynamics equations are reversible in time. To achieve propulsion, 
microswimmers must therefore deform in a way that is not invariant under time reversal. Here, 
we investigate dispersal properties of the micro-alga Chlamydomonas Reinhardtii, by means of 
microscopy and cell tracking. We show that tracked trajectories are well modeled by a correlated 
random walk. This process is based on short time correlations in the direction of movement called 
persistence. At longer times, correlation is lost and a standard random walk characterizes the trajec- 
tories. Moreover, high speed imaging enables us to show how the back-and-forth motion of flagella 
at very short times affects the statistical description of the dynamics. Finally we show how drag 
forces modify the characteristics of this particular random walk. 

PACS numbers: 47.63.-b, 47.57.-s, 47.50.-d 



Cell motility [1 is crucial to many biological pro- 
cesses including reproduction, embryogenesis, infection, 
etc. Many microorganisms are able to propel themselves, 
bacteria, sperm cells, microalgae, etc. A quantitative un- 
derstanding of the hydrodynamics of flagella and cilia is 
thus of great interest |2j [3] . 

One of the peculiarities of the swimming of microor- 
ganisms is that it occurs at very low Reynolds num- 
bers which is very different from our usual experience 
of swimming at our meter length scale [H [5]. Indeed 
when inertia is negligible as compared to viscous forces 
(i.e. Reynolds number Re is lower than unity), in order 
to achieve propulsion, swimmers must deform in a way 
that is not invariant under time reversal. This is known 
as Purcell's scallop theorem [4]. In living systems, sev- 
eral different strategies are used to achieve propulsion in 
such conditions: the E. coli bacterium uses a rotating 
flagellum at the "back" of its body, sperm cell propul- 
sion relies on the asymmetry of their flagellar bending 
waves, the power and recovery strokes of the two front 
flagella of Chlamydomonas Reinhardtii (CR) algae are 
asymmetrical. 

Flagellar propulsion in CR induces complex swimming 
behavior of cells. Over short time scales, the cells un- 
dergo an oscillating movement with changes in velocity 
direction occurring at the same frequency as the beating 
frequency of flagella. On a time scale longer than the 
period of beating, average swimming behavior is direc- 
tional. Eventually, on larger time scales, direction is lost 
and swimming trajectories resemble a random walk. 

CR is a 10/im motile bi-flagellated unicellular alga. 
The cell is spheroidal in shape with two anterior flag- 
ella [6]. It belongs to the puller type of swimmers as it 
uses its front flagella to propel itself, producing what a 
breast stroke-like movement. The swimming direction of 
the cells can be controlled by stimulus gradients, a phe- 
nomenon known as taxis, such as chemotaxis, rheotaxis 
or phototaxis. Gradients are not used in our experiments 
in order to avoid any external tropism on the motility. 
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FIG. 1: (color online) (a) Probability distribution functions 

of angles = arccos[/c(£o)-&(£o + ^)] ^ or different times t rang- 
ing from 0.05£ c to 3.33£ c . Here £ c =2.3s and the viscosity of 
the medium is 2.4mPa s. Only few data points are displayed 
for the sake of clarity, (b) Correlation functions of direction 
C(t/t c ) as defined in the text (Log-Lin scale). Time has been 
rescaled by the decaying time of an exponential. The differ- 
ent symbols correspond to the viscosities used by varying the 
concentration of dextran. 



Wild-type strains were obtained from the IBPC lab in 
Paris [7 J. Synchronous cultures of CR were grown in 
a Tris- Acetate Phosphate medium (TAP) using a 12/12 
hour light /dark cycle at 22° C. Cultures were typically 
grown for two days under fluorescent lighting before the 
cells were harvested for experiments. 

We studied the swimming dynamics of this microor- 
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ganism by means of bright field microscopy imaging on 
an Olympus inverted microscope coupled either to a CCD 
camera (Sensicam, Photon Lines) used at a frame rate of 
10Hz or to a high speed CCD camera (Miro, Phantom) 
used at a frame rate of 400Hz. Long time experiments 
used a xlO magnification lens whereas we used a x64 
lens for high speed imaging. Two hundred micrometer 
thick glass chambers were coated with bovine serum al- 
bumine to prevent cell adhesion. The imaged cells were 
located 30 to 60 micrometers from the glass walls. A 
red light filter was used in order to prevent phot ot axis. 
Cell tracking [8 was performed using IDL (Interactive 
Data Language) with a submicron precision in the detec- 
tion of hundreds of cells (high speed experiments) and 
thousands of cells for long time sequences. To quantify 
the effect of drag on the cell dynamics, small amounts of 
short chain dextran (Sigma Aldrich) were added to the 
culture medium. The chains were short enough for non- 
Newtonian effects to be absent and long enough to avoid 
damaging the cells with osmotic effects. This allowed 
the viscosity r] of the medium to be varied between 1.5 
and 3.7mPa s. The range of viscosity is restricted to this 
interval to ensure the viability of the cells. 

Let us first recall the global dynamics of swimming, 
i.e. over time scales of the order of a few seconds. Cell 
trajectories are found to be correctly modelled by a per- 
sistent random walk [9-11 . Cells swim in an almost fixed 
direction for a typical time of about one second. This 
stage corresponds to a ballistic regime characterized by 
a mean velocity V. The ballistic regime ends when the 
swimmers make a turn. A new direction of motion is then 
observed due to the desynchronization of the pair of flag- 
ella [12]. At long time scales, the dispersal properties of 
the swimmers are random-like [I3j [14] . To describe this 
specific random walk quantitatively, we measured differ- 
ent statistical quantities of interest. Let us first define 
a persistence angle 6{t) = arccos[fc(to)-fc(^o where 
k(t) is a unitary vector in the direction of movement 
at time t. Thus, a value of close to zero reflects a 
certain persistence of the trajectory. We measured the 
probability distribution function of angles for different 
times t. At short time scales, angle distribution peaks at 
around zero, characterizing the directional persistence in 
swimming trajectories. Over longer times, we observed a 
broadening of the distribution and eventually we ended 
up with an equidistribution of angle values characteristic 
of a random walk (figure [TJa). This phenomenon is even 
better quantified by looking at the mean value of angle 
distribution or equivalently at the correlation function of 
direction defined as C(t) =< k(to).k(to + t) >, where 
<> is an average over time to and over all tracked tra- 
jectories. Correlations with infinite decay time (C(t) = 1 
for all t > 0) correspond to direction correlations pre- 
served over arbitrarily long times, i.e. a purely ballis- 
tic regime; whereas a zero life-time (C(t) = for all 
t > 0) corresponds to standard random walk behav- 
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FIG. 2: (color online) Mean Square Displacements < r 2 (t) > 
of cells as function of time for different viscosities of the 
medium. The different symbols represent different bath vis- 
cosities. The legend for symbols is the same as in figure [3^b). 
The solid lines represent a slope 2 and the dotted line a slope 
1 in a log-log scale. Typical 2D trajectories of a few cells are 
displayed in the insets. In the top left inset, cells are swim- 
ming in the nutritive medium (viscosity r/ = 1.5mPa s), in the 
lower right corner, the medium is rich in dextran (77 = 3.7mPa 
s). Trajectories are represented on the same scale for better 
comparison and both lasted 0.5 s, their starting positions were 
all shifted to the origin. 

ior (figure [TJb). The correlation functions decay expo- 
nentially over a characteristic time t c . This correlation 
time t c is related to the mean time of persistence over 
which the direction of swimming is preserved. The dif- 
ferent symbols in figure [l]b correspond to experiments 
where the concentration of short chain dextran was var- 
ied, hence modifying the viscosity of the medium from 
1.5mPa s to 3.7mPa s. As viscosity increases, correla- 
tion time t c increases (data not shown) from 1.5 to 3.9 
seconds. 

The global dynamics of swimming of CR can thus be 
described as a correlated random walk characterized by 
a ballistic regime (with a mean velocity V) and a decor- 
relation process (over a characteristic time t c ) due to the 
turns made by the cells. As a consequence, a persistence 
length C is naturally defined as the product Vt c . From a 
statistical point of view, such a bahaviour is described by 
the mean square displacement of cells < r 2 (t) > which 
is linear for long times (t^> t c ) and quadratic at shorter 
times (t < t c ) p~3j[l4]. At even shorter times, the dynam- 
ics reflect the consequences of low Reynolds swimming, 
i.e. a non-reciprocal movement of flagella. This then 
leads to a zigzagging motion of cells due to the back- 
and- forth motion of flagella [T5] . 

In the present case, cells of diameter 2R ~ 10/im 
are moving at a velocity V around 50/im/s in a wa- 
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FIG. 3: (color online) (a) Probability distribution functions 
of angle at short time scales ranging from t = 1/ (28/) tot = 
1/ f. In these experiments, cells are swimming in a medium 
of viscosity r/ = 2.4mPa s so that their beating frequency is 
/ = 14.3Hz. (b) Correlation functions of direction as defined 
in the text. Time has been rescaled with the period 1/ f of the 
signal which corresponds to the beating frequency of flagella. 
Symbols correspond to different viscosities of the medium. 
The solid line represents the function cos(27r£#) exp(— £#). 



ter like medium (viscosity rj ~ ImPa s and density 
p ~ 10 3 kg/m 3 ). This represents a very low Reynolds 
number of the order of Re = pVR/r] ~ 2.5 x 10 -4 . The 
propulsion strategy of CR consists in swimming in a kind 
of breaststroke where the pair of flagella are wide open 
during the forward movement and folded along the cell 
body during the backward movement [16 . Hence, vis- 
cous friction is high when the pair of flagella are fully ex- 
tended during the forward movement and friction is lower 
during the backward movement. The symmetry under 
time reversal is thus broken and propulsion is ensured. 
However, because inertia has no role in this regime, this 
kind of propulsion leads to a back-and-forth movement 
of the cell in which the velocity is alternatively positive 
and negative. High speed imaging (400Hz) allows us to 
resolve the very short time dynamics due to flagella beat- 
ing and thus to study the consequences of this back-and- 
forth movement on the properties of the swimmers' ran- 
dom walk. 

The insets in figure [2] show typical cell trajectories im- 
aged at 400Hz: the back-and-forth movement of swim- 
mers due to the absence of inertia (Re << 1) together 
with the long time swimming behavior are visible. In 



these examples, the cells are swimming either in a nutri- 
tive medium of viscosity r] = ImPa s (top left inset) or 
in a dextran-rich medium of viscosity 3.7mPa s (bottom 
right inset). Cells have a net forward movement corre- 
sponding to the power stroke, followed by the recovery 
stroke that propels the cell backward. As the distance 
traveled forward is longer than the backward movement, 
the cells ultimately progress forward. However, these 
fluctuations in the direction of the velocity have conse- 
quences on the measured statistical quantities [17 that 
we will discuss now. 

The measured mean square displacement < r 2 (t) > 
shows a plateau region at very short time (t <C t c ) that 
reflects the transition between two quadratic regimes : 
on the one hand, a fast ballistic regime characterized by 
the instantaneous velocity u of swimmers and, on the 
other hand, a slower ballistic regime corresponding to 
the mean velocity V of swimming which is the resulting 
forward velocity over several back-and-forth movements. 
The position of the plateau therefore corresponds to the 
beating frequency / of the swimmer, which depends on 
the viscosity of the surrounding medium. To quantify the 
back-and-forth swimming motion of the cells, we mea- 
sured the angle probability distribution function. Fig- 
ure Oa shows distribution functions for different times. 
For a given short time t, the distribution of angles 6(t) as 
defined earlier, peaks at around zero, reflecting a given 
direction at very short time. For longer time scales (close 
to l/2f), anti-correlation in cell direction resulted in new 
distribution peaks at values of around ±7r. When a new 
stroke is produced, the measured angle is again close to 
zero giving rise to a peak around zero. Hence, angle dis- 
tributions have a periodicity which reflects the beating 
frequency. This is shown in figure [3ja as the distribu- 
tions are very similar at times shifted by 1/(2/), where 
the typical frequency of the beating / is deduced from 
the periodical nature of the correlation function of di- 
rection. Figure [3]b shows such correlation functions at 
varying /xt, the product of time multiplied by the fit- 
ted frequency of the signal. Data are well described by an 
exponentially attenuated cosine function. The different 
symbols correspond to different viscosities of the medium. 
The exponential decay of the correlation function should 
reflect the turns in direction the cells eventually perform 
within a characteristic time t c . However, due to the 3D 
nature of the trajectories and the 2D geometry of our 
setup, correlation was attenuated faster than that. 

The other consequence of the fact that swimming is 
produced at low Reynolds number is that propulsion re- 
quires non-zero drag forces. Viscous friction is thus cru- 
cial in the dynamics of micros wimmers. By varying the 
viscosity of the medium, we were able to draw some con- 
clusions about the effects of friction forces on the loco- 
motion of microorganisms such as this micro-alga. 

Here, short time dynamics of swimming can be fully 
described by few mean quantities: flagella frequency 
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FIG. 4: (a) Beating frequency, obtained from the period of 
oscillations in the correlation functions of directions, versus 
the inverse viscosity of the medium, (b) mean modulus of 
velocity u as a function of the inverse viscosity of the medium. 
Solid lines represent linear regression giving: r\f — 0.045 =b 
O.OlPa and rju = 0.15 ± 0.04Pa./xm 

beating / (deduced from a cosine fit in figure [3]b) and 
the mean modulus of instantaneous velocity u, which is 
the velocity achieved during a power or a recovery stroke. 
We studied the effects of viscous forces on these quanti- 
ties. Velocities and beating frequency are found to be 
inversely proportional to the viscosity of the bath (fig- 
ure [4]). As viscosity increases, the beating frequency 
decreases, varying from 30Hz to 13 Hz with a viscos- 
ity variation from 1.5 to 3.7mPa s (figure [4^ a) giving a 
slope rjf = 0.045 ± O.OlPa. Accordingly, velocity de- 
creases from 135 to 75/im/s (figure |4^b) giving a slope 
rju = 0.15±0.04Pa /im. These results supports the idea of 
imposed- force locomotion [14 j. The corresponding stall 
force, which is proportional to the product 77 x u, is then 
constant. 

The velocity u can be related to the mean propulsion 
force on the cell body by Stokes' law. Let's now assume 
that a power stroke (respectively a recovery stroke) re- 
sults from the friction length £_l (resp. £y) of the flagella 
moving perpendicular (resp. parallel) to its long axis. 
Moreover, the beating frequency / can be related to the 
friction of flagella acting on a typical distance of one 
cell diameter: 6ttt]Ru = 2Rf r q(£ ) ± + £n). Using measure- 
ments of the slopes rju and 77/, we can estimate a sum 
£_l + £|| ~ 32/ira. Using the friction coefficient expres- 
sions of a cylinder given in [18], this leads to an aspect 



ratio of 200 for a 10/xra long flagellum with a radius of 
25nm. This is a reasonable estimate [6 considering that 
flagella are not exactly perpendicular and parallel to the 
flow during power and recovery strokes. 

In this work, we quantified the complex dynamics 
of swimming at short time scale using high speed mi- 
croscopy imaging and particle tracking techniques. This 
study allowed us to analyze the breaststroke like swim- 
ming of a CR cell in a fluid at low Reynolds number 
and how this swimming is influenced by the viscosity of 
the ambient fluid. It showed how the friction acting on 
a CR cell can be qualitatively extracted from the back- 
and-forth motion of a thin and elongated pair of flagella. 
This means that a description in terms of time averaged 
flows is not to be encouraged for such systems [3j [19] . 
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